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ABSTRACT 

Over the past decade, many works on the modeling of wire¬ 
less networks using stochastic geometry have been proposed. 
Results about probability of coverage, throughput or mean 
interference, have been provided for a wide variety of net¬ 
works (cellular, ad-hoc, cognitive, sensors, etc). These re¬ 
sults notably allow to tune network protocol parameters. 
Nevertheless, in their vast majority, these works assume that 
the wireless network deployment is flat: nodes are placed on 
the Euclidean plane. However, this assumption is disproved 
in dense urban environments where many nodes are deployed 
in high buildings. In this paper, we derive the exact form of 
the probability of coverage for the cases where the interfer- 
ers form a 3D Poisson Point Process (PPP) and an approx¬ 
imation for the 3D Modified Matern Process (MMP). We 
compare the 3D model with the 2D model and with sim¬ 
ulation results. We comment the adequacy of each model 
in function of the parameters of the nodes (emission power, 
reception threshold, MAC protocol, etc.) and the height of 
the buildings in the simulations. 

Keywords 

Urban Wireless Networks; 3D; Stochastic Geometry; CSMA 

1. INTRODUCTION 

Stochastic geometry have been largely used to study and 
design wireless networks, because in such networks the in¬ 
terference and thus the capacity is highly dependent on the 
positions of the nodes . Stochastic geometry indeed 

allows to take into account the spatial component for the 
analysis of wireless systems performance at a very low com¬ 
putational cost (in several cases closed form expressions are 
available [^). Nevertheless, most of the works in the liter¬ 
ature focus on networks deployed on the Euclidean plane. 
The motivation for this work comes from the intuition that 
for dense urban networks such as WiFi private networks, 
sensor networks, connected objects deployed in dense ur¬ 
ban areas, modeling the network by projecting all the nodes 


PREPRINT SUBMITTED TO SPRINGER . 


on a plane will lead to an inaccurate representation of the 
interference suffered by a node. We think that with the ex¬ 
plosion of the number of devices expected in the ISM bands, 
notably because of smartphones, tablets, connected objects 
(Internet of Things) and sensors, there is a need for an accu¬ 
rate modeling of the interference suffered by such devices in 
dense urban areas. Indeed, accurate interference modeling is 
essential in order to design efficient mitigation mechanisms 
at different layers (medium access schemes, coding schemes, 
retransmission mechanisms, robust routing, etc.). In order 
to be accurate, the third space dimension has to be taken 
into account. 

Throughout this paper, we model the node positions as 
a Poisson Point Process (PPP). We argue that the PPP 
seems to be a good model for dense urban areas where many 
networks are concurrently deployed in the same frequency 
bands (notably ISM bands) in a chaotic manner. Unlike 
cellular networks, for WiFi, sensor and loT networks no 
planned deployment can be assumed in general [^. 

For the numerical applications, we take the example of 
private WiFi networks in a dense urban area such as the 
center of Paris. We consider, the V7e arrondissement of 
Paris which contains around 32500 dwellings and covers 
a surface of 2.15km^. We assume that each dwelling has one 
WiFi router. In the 2D case, they form a Poisson Point Pro¬ 
cess (PPP) on the plane. In the 3D case, they are also spread 
on the z axis. If we assume a mean building height of 20 me¬ 
ters for the 3D case, it gives intensities of A = 1.51 x 10“^ 
nodes per and p = 7.56 x 10“^ nodes per respec¬ 
tively for the 2D and 3D processes. These intensities may 
seem high, but they actually do not take into account the 
devices (laptops and smartphones for instance) connected to 
the routers which also increase the interference. 

In the remainder of this paper, we derive the probability 
of coverage, first with no medium access control (the results 
are applicable to ALOHA [^), and then with a CSMA access 
when nodes are distributed according to a 3D PPP. 

The main contributions of this paper are as follows: 

• we give an approximation of the probability of coverage 
for CSMA access in the 3D case; 

• we compare 2D, 3D analytical models and simulation 
in order to determine the most suited model in function 
of the parameters of nodes and network topology (no¬ 
tably the height of the deployment along the z axis). 

The paper is organized as follows. Section presents and 
comments related works on stochastic geometry and 3D net¬ 
work modeling. In Sectionj^ we derive the exact form of the 


probability of coverage for PPP interference for arbitrary 
emitter-receiver distance and compare the results with the 
2D model and simulations. It gives us insights on the prob¬ 
lem for a simple case. In Section]^ we present the Modified 
Matern Process used to model CSMA access and derive an 
approximation of the probability of coverage for the 3D case. 
We then compare this approximation with the one of the 2D 
case and simulations. Section[5]concludes the discussion and 
provides planned future developments. 


2. RELATED WORKS 

Even if most of the research efforts are focused on 2D 
networks, 3D networks have been emerging during the past 
few years. It is notably the case for WSNs |22[ under¬ 
water mobile networks Unmanned Aerial Vehicle (UAV) 
networks [21| . 

From a more theoretical point of view, 3D networks have 
been investigated in terms of capacity [12[ |17| and scaling 
laws have been provided. Nevertheless, in the present work 
we focus on the stochastic geometric approach to the study 
of wireless networks in the sense of [7 14 . This approach 
has the advantage to provide tractable results for the prob¬ 
ability of coverage in the case of the PPP model and good 
approximations for other models [^. 

Whereas, many theoretical works focusing on stochastic 
geometry for wireless networks consider dimension d , 

when it comes to applications to specific cases, the chosen 
space is nearly always the Euclidean plane [18[ . The 

only work explicitly covering the 3D case, to the best of 
our knowledge, is the recent [^. In this work, the authors 
derive the probability of coverage with the two slopes prop¬ 
agation model for a 3D PPP. This work considers cellular 
networks, the receiver is thus attached to the nearest point 
of the PPP (the base station with the highest average re¬ 
ceived power). In our work, we first focus on a 3D PPP 
model as well, but with a single slope propagation model 
and, more importantly, the receiver is receiving from an ar¬ 
bitrarily placed emitter. We argue that this model is more 
relevant in non operated ISM band networks such as WiFi, 
loT, WSNs, etc, because the receiver cannot always connect 
to the closest node. Indeed, the closest node might not be¬ 
long to the same network. Moreover, in our work, we mainly 
use the simple PPP model to easily compare simulation and 
theoretical models for different building heights (maximum 
value on the 2 axis) and thus evaluating the relevance of the 
3D approach. We also then move to the Modified Matern 
Process which is a more realistic model CSMA networks. 

Many works modeling CSMA access through stochastic 
geometry have emerged in the literature these past ten years 
|18[|16[[l0| . These works are based on a modihed version of 
the Matern Point Process for which the points cannot live 
too close to each others. This allows to model the contention 
radius of CSMA as follows. The nodes of the network form 
a PPP on the Euclidean plane and they contend to access 
to the medium. Each node picks uniformly a value between 
zero and one, the node is selected in the process if it does 
not detect any node with a smaller mark value (this models 
the backoff procedure). The resulting process is the Modi¬ 
fied Matern Process (MMP) which will be detailed in section 
m As there is no known exact formulation of the interfer¬ 
ence in such process, [18[ |16| and ]10| use approximation 
techniques. They show that their approximations lie close 
to simulations results. We extend these works by consider¬ 


ing 3D distribution of the nodes and compare the resulting 
model to simulation. To the best of our knowledge, this has 
not been considered in the literature. 

In the remainder of this paper we show how the proba¬ 
bility of coverage is changed when going from 2D to 3D for 
Poisson distributed nodes in the simple case where all nodes 
can be interferers, and in the case a CSMA access protocol 
is used. We compare the theoretical prediction with simula¬ 
tions. We also show that going from 2D to 3D is not trivial 
especially in the CSMA case because changing the dimen¬ 
sion dramatically affects the form and thus the tractability 
of the expressions. 


3. 3D POISSON POINT PROCESS INTER¬ 
FERENCE 

In this section we derive the probability of coverage for a 
terminal at distance d from the emitter when the interferers 
are distributed according to a PPP 4? in R®. Whereas it is 
already given for R" and for R^ in [^, here we find inter¬ 
esting to emphasis the differences between the 3D and 2D 
cases and to show that considering the latter can lead to 
inaccurate representation of the interference. We consider 
an interference limited network. We thus use the Signal to 
Interference Ratio (SIR), which is dehned as follows: 


SIR = 


hd- 


( 1 ) 


with h the fading coefficient between the emitter and re¬ 
ceiver, d the emitter-receiver distance, a the pathloss expo¬ 
nent and I = interference where Qi is the 

fading coefhcient between the interferer Xi € "F and the re¬ 
ceiver and Vi the distance between them. As in 15 


assume that all the nodes emit with the same power so it is 
simplihed in the expression of the SIR. We have to note that 
unlike many works in the literature, we do not consider that 
the receiver is connected to the nearest node of the PPP 
because it is not always the case for the type of networks we 
consider (private WiFi is a good example) so the probability 
of coverage depends on the emitter-receiver distance. The 
probability of coverage is defined as the probability that the 
SIR is over a given threshold /3: 


P,{p,P,a,d) = P{SIR>li} 


( 2 ) 


Its expression, when considering Rayleigh fading between 
the emitter and receiver, can be derived by the classic ar¬ 
gument notably found in and is equal to the Laplace 
transform of the interference shot noise in interference lim¬ 
ited networks: 


P,{p,P,a,d) = p[^^> 

= Ci{pPd’^) (3) 

(a) follows from H ~ exp{p). 

The expression of the Laplace transform for the 3D PPP 
case is as follows: 
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(b) follows from Gi being i.i.d. and independence with 
-Ri, (c) fol low s from the probability generating functional of 
the PPP [20| , (d) follows from the MGF of Gi ~ exp{fi) 

and (e) with the change of variable u = ^ ^ differs 

-A7rd2^2/“ ^/°( ^^)dv 

from the 2D case which is given by e o ' 

We note that the solution of the integral in the general case 
(dimension D) is: 


+ 00 
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( 5 ) 


with 2 F 1 the Gaussian hypergeometric function. This 
function gives highly different values even for the same pathloss 
exponent a in the 2D and the 3D cases. The other main dif¬ 
ference is the exponential decay which is in dP in the 2D case 
and dP in 3D. We give the closed form for a = 4 for the 3D 
case: 


Fc(p,/3,4,d) 


g-p(4/23/2).^2jj3^3/4 


( 6 ) 


emitter. Given the steepness of the decay, the error might 
be important: for instance, there is 0.23 absolute difference 
at 2 meters. We note that the Pc is getting very rapidly 
low (after few meters Pc goes to zero) which is not realistic. 
We observe this because, in this simple case, all the nodes 
transmit at the same time, whereas it is not the case in re¬ 
ality where a MAC protocol is used. In Section]^ we focus 
on the modeling of a more realistic setup: a 3D CSMA net¬ 
work. Nevertheless, We can note that expression § is valid 
for modeling an ALOHA medium access. In this case, the 
intensity of the process should be multiplied by the proba¬ 
bility for a node to emit . ALOHA access is currently one 
of the considered access mechanism in the loT standard Lo- 
RaWAN from the LoRa Alliance . The 3D PPP model 
could thus be used to model such loT networks in dense 
urban areas. In the next subsection, we present simulation 
results which advocate for the relevancy of the 3D model. 

3.1 Comparison with simulation results 

In this section, we compare the theoretical and simulation 
results. We evaluate the relevancy of the 3D approach for 
dense urban areas by considering the following question: for 
which building height does it makes sense to model the in¬ 
terference process as a 3D PPP? Indeed, in reality, the nodes 
will be spread within a finite height because of the limited 
height of the buildings, so a border effect may appear com¬ 
pared to the theoretical model. 

The simulation setup is as follows: the point process is 
generated in box of 200mX 200mX Z with Z € 10m, 50m, 100m 
The 3D intensity of the process is calculated based on the 
2D intensity: p=F reason for this is that in the 2D 

case we project all the points on a plane, the intensity is thus 
maximal. In the 3D case the points are also spread along 
the 2 ; axis. Then for each realization of the PPP, the SIR 
is computed at the origin (the center of the aforementioned 
box) thanks to equation 0- The simulation and model pa¬ 
rameters are the following: A = 1.51 x 10“^ (we keep the 
values from the example of Section]^, /3 = 10 and a = 4. 
The emitter-receiver distance d varies from 0 to 10 meters. 



Figure 1: Probability of coverage for 2D and 3D PPP inter¬ 
ference 

In Fig. we plot Pc{p, P,4,d) for the 2D and 3D cases 
with A = 1.51 X 10“^ and p = 7.56 x 10“'* (as motivated in 
Section]^ and /3 = 10, the plot is in function of the emitter- 
receiver distance d. We observe that Pc, in the 2D case, 
is highly underestimated when the receiver is close to the 



Figure 2: Comparison of theoretical and simulation results 
for PPP interference 

Fig. [^depicts the comparison of the probability of cover¬ 
age Pc between simulation and theoretical results for dif¬ 
ferent building heights (parameter Z) in function of the 
emitter-receiver distance. For each simulation point, 10 000 


















PPP realizations are produced. The 95% confidence inter¬ 
val is plotted, but is barely visible (it is very small). The 
2D model corresponds to the solid curve, we observe that it 
underestimate the probability of coverage compared to 3D 
models and simulations. The three dashed and dotted curves 
correspond to the 3D theoretical models. They are obtained 
from expression ^ with p = ^ and Z € 10m, 50m, 100m. 
We observe, as expected, that the probability of coverage is 
higher when the considered height Z is greater because the 
nodes are more spread and thus the distances to the inter- 
ferers are larger on average. For Z = 50m and Z = 100m 
the simulations match the theoretical prediction very well. 
Nevertheless a tiny deviation can be noticed at the tail of 
the curves. This deviation is larger for Z = 50m than for 
Z = 100m. In the Z = 10m case, the deviation of the the¬ 
oretical model from the simulations is very significant. We 
can interpret these deviations as a border effect when the 
height of the box is reduced: the 3D theoretical model as¬ 
sumes no limits on the x, y and a axis. We also observe 
that even when the height of the box is reduced to Z — 10m 
the 2D theoretical model is not a good representation of the 
simulation data either. In this case, we note that the 3D 
model is a better match for the start of the curve and the 
2D model is better for the tail of the curve. 

We conclude that the 3D model seems to be accurate 
enough when the box height is sufficiently large to avoid 
border effect. Moreover, even when the height is small (the 
Z = 10m case) 3D model still provides a good representa¬ 
tion of the probability of coverage close to the emitter be¬ 
cause the border effect is larger for longer emitter-receiver 
distances. 

In the next part of this paper, we consider another medium 
access model: CSMA. More precisely, we extend to R® an 
approximation of the probability of coverage for the MMP 
model. We then compare the theoretical model predictions 
to simulation results. 


4. 3D MODIFIED MATERN PROCESS IN¬ 
TERFERENCE 

In this section, we treat the 3D case for the Modified 
Matern type II Process (MMP) which is considered in many 
works [18[ |10| in order to model CSMA access. With this 
model, nodes too close to each others cannot transmit at the 
same time, interference is thus reduced. The MMP model is 
detailed in Section [4. II 

Since the exact Laplace functional is not known for the 
MMP, the technique used in the previous section for PPP 
(notably step (c) in equation 0). cannot be applied di¬ 
rectly. In previous works, several different approximation 

1 10]. In this paper, we use 


techniques are used 18 


the technique which consists in: 

1. assuming a finite region containing the potential con¬ 
tenders for a node; 

2. deriving the resulting MMP intensity pcsma', 

3. considering only dominant interferers; 

4. approximating the MMP with a PPP of intensity pcsma 
outside of the emitter contention domain. 

This approximation technique is notably used in [10| for 
an alternative version MMP. In our work, we derive expres¬ 
sions for the 2D and 3D cases, whereas jlOj is only for the 


2D case and the alternative MMP. This technique allows 
to obtain an expression of the probability of coverage for 
the MMP, despite the fact that the exact expression is not 
known. It is detailed for the 3D case in the remainder of 
this section. 

4.1 The Modified Matern type II Process 

In this part of the paper, we use a modified version of 
the Matern type II process [^. The classic Matern type II 
process is built from a marked PPP: each point of the PPP 
is marked with a real m € [0,1] and the PPP is thinned by 
retaining only the points which have the smallest mark in 
a ball of radius D centered on themselves. Formally, 4?^ = 
{xi £ <l>|mi < mjt/xj £ ^r\B{xi,D)\xi], with the PPP 
in R®, rrii the mark of point Xi and B{xi, D) a ball of radius 
D centered on Xi. In the modified version of the process, 
there is no fixed radius D, but the radius is replaced by the 
notion of detection: a node from the marked PPP is kept in 
the MMP if it does not detect any signal from a node with 
a lower mark. The signal from a node is detected if it is 
above a threshold Td. So formally the process is defined as 
■hm = {xi £ 43|mi < mj,\/xj : Pthijr~Z > Ta}, with Pt the 
transmission power, hij and Vij the fading coefficient and 
distance between Xi and Xj respectively. 

The MMP allows to model the position of the nodes ac¬ 
cessing the medium concurrently with a CSMA MAC scheme. 
The underlying PPP represents the positions of the nodes. 
In the following subsections, we detail the steps of the afore¬ 
mentioned approximation technique applied for the 3D case. 

4.2 Defining a finite region for contenders 

In reality, the region which contains the contenders of a 
given node (the nodes which can be detected by that node) 
presents an irregular shape because of the randomness of 
the fading coefficients. Nevertheless, as mentioned previ¬ 
ously, the first step of the approximation technique consists 
in considering a fixed radius, called the detection radius (de¬ 
noted Xd), outside of which the probability (denoted ed) that 
a node is detected is arbitrary low. Thus the contenders of 
a node Xi are contained in B(xi,rd) which is a ball in the 
3D case. It allows to simplify the expressions without losing 
too much accuracy [10| . Formally, the detection radius is 
defined as follows: 

P{Pthijr'^°‘ > Td) <ed<^rd = 

^rd = 

if we consider Hij ~ exp{p). 

4.3 The intensity of the Modified Matern Pro¬ 
cess 

The intensity of the resulting process (considering only 
contenders in B{xi,rd)) is given as pcsma = pPcsma [20] , 
with Pcsma the probability for a node of the underlying PPP 
to be retained and p the intensity of the PPP. Thus Pcsma is 
the probability that among the contenders of a node Xi (the 
nodes of <I> which are inside B{xi, Vd) and that are detected 
by Xi) none have a smaller mark than Xi. 

Formally, it is given as: 
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( 8 ) 


The details of the argument are similar as those of or 
[10| and thus they are not reproduced here. P„ and Pd are 
respectively the probability to have n nodes in B{xi, r^) and 
the probability to detect a node which is in B(xi,rd)- Pn is 
thus the probability that there are n nodes in a volume of 
(4/3)7rrj^ which is given by the Poisson distribution. For Pd, 
in the 3D case we have: 


Pd = PiPthijT-;^ > Td) 

= Ep.^- 




/ 


3r: 


3 


Pt 




3 /q 


( 9 ) 


with r(a) and r(a,b) the Gamma function and incomplete 
upper Gamma function, (f) follows from Pij ~ exp{p,) and 
fRijiUj) is given by: 


fRij 


Of 


( 10 ) 


fRij {fij) is the density function of the random distance 
between Xi and Xj in B{xi,rd)- In order to prove that ex¬ 
pression m is correct for the 3D case (for the 2D case it 
is given in |19|) we have to remember that Xj is uniformly 
placed in B{xi, Xd) because of the underlying PPP proper¬ 


ties, so we have FR,Arij) = 


{4/3)irr3 ' 


Taking the derivative 


with respect to Xij yields the result. Expression § differs 
from the 2D case where fRi-ifij) is 19 , but the ex¬ 


pression remains reasonably easy to compute in the 3D case 
(the Gamma functions can be efficiently computed). 


From ^ and Pcsma = pPa 
tensity of the MMP is: 


we conclude that the in- 


I _ g-p(4/3)7rr|p,j 

= {4/3)7vrlPd 

4.4 Dominant interferers vulnerability radius 

Once we have obtained the intensity of the MMP, we con¬ 
sider only the interferers which can corrupt the signal on 
their own and we assume that they all lie in a region around 
the receiver, called the vulnerability region. Similarly to the 
detection radius, we define the vulnerability radius as the 
radius for which nodes beyond that limit have an arbitrary 
low probability (again noted e„) to make the SIR at the re¬ 
ceiver drop under the reception threshold. The receiver is 
placed at the origin and noted o. The test emitter is part of 
the process and noted Xi. Formally the vulnerability radius 
is defined as follows: 
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( 12 ) 


with - 77 ^ the ratio of exponential random variables of pa- 

PI jo 

rameter p so Fhj^ (1) = 1 — . 


4.5 Probability of coverage 

We give an approximation of the probability of coverage 
^pcsma^ for a typical receiver at the origin. As previously 
mentioned, the MMP is approximated with a PPP of the 
same intensity outside of the contention domain of the emit¬ 
ter Xi, this approximation is known to be very close to simu¬ 
lation results [^. In this context, the probability of outage 
is the probability that among the nodes in the vulnerability 
radius which coexist with Xi (outside of the contention re¬ 
gion of Xi) some will be able to make the SIR drop under the 
threshold /3. The probability of coverage is the complement 
of the probability of outage: 
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with Pk 




K,. 


i(4/3)7rrj, Pd' 


and Pp are respectively the probability that a node coexists 
with the emitter and the probability that the SIR drops 
under the reception threshold. Pd' and Pp are derived below 
in this section. The full detail argument for a result close to 
Equation (131 (2D case) can be found in (equation (131 
derivation only involves well-known convergent series). Deep 
differences from the 2D case are found in the expressions of 
Pd' and Pp derived below (and also in expression (111 of 
pcsma as shown in Section [T^ . 

Pd' is the probability for a node in S(o, r„) to be in the 
contention domain of Xi (to be detected by Xi): 
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where 2^1 is the Gaussian hyper-geometric function, (g) fol- 

3r^ 

lows from F(1) = 1 — 3 ^ and fRj^{rjo) = and (h) 


from the change of variable w = 
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c = 


r(2/a,B), D = r(3/Q,B), E = r(4/a,_B), F = r(2/Q,A), 
G = r(3/a, A) and H = r(4/a, A). 


Here the expression of /_R; . (r^j) changes from (10 1 in Sec¬ 


tion 4.3 because we are no more interested in B[xi, Vd) but 
rather in B{o,ry). As can be seen in Fig. [^we now have 
two cases: either B{xi,rij), the dotted sphere, is contained 
in B{o,rv), the dashed sphere (0 < rij < r„ — Vio), or it is 
not (r„ — Vio < Tij < + rio). In the first case, the expres¬ 

sion is similar to ( |10[ | . In the second case, the expression for 
fRij (rij) changes because we know that Xj cannot be out¬ 
side of B{o, Tv) so it must lie in the intersection of B{xi, Vij) 
and B{o, Tv)' 


fRiiirij) 


' 3 r • - 

^, for 0 < Tij < r„ - Vi, 
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for r„ — Vio < n 


(15) 

Interestingly, we note that the second piece of the expres¬ 
sion of /r.^ ifij) is rational whereas in the 2D case it includes 
trigonometric functions (arcsine) 10 . It allows to have an 
exact expression for Pd' in the 3D case (expression (Et) 
whereas, to the best of our knowledge, numerical integra¬ 
tion has to be used in the 2D case. 

Pp is the probability that a node in B{o, r„) is able to make 
the SIR at the receiver drop under the reception threshold 

P-- 


Figure 4: Probability of coverage for 2D and 3D MMP 

Fig. [^depicts the probability of coverage for the 2D and 
3D cases for the MMP. (Equation ( |13[ )) is given in 

function of r^o, the emitter-receiver distance. The curves 
are plotted with the values Pt = lOOmW, Td = —76dBm, 
/3 = 10, a = 4, /r = 1, Cd = 10~®, e„ = 10“^, and p and 
A from Section First we have to note that in this case, 
the communication range is more realistic (even if in real- 
+ Tty in dense urban environment, we may encounter pathloss 
exponents greater than 4). We observe that, with the afore¬ 
mentioned parameter values, the probability of coverage is 
strictly higher in the 2D case for emitter-receiver distances 
above 50m. We can interpret this result as follows: for 
the PPP in Section]^ we remarked that the 3D probability 
of coverage tends to drop later than the 2D, but then the 
slope is steeper. For the MMP model, the difference in the 
drop position is suppressed by the contention mechanism 
(it avoids close interferers), but the difference in the slopes 
remains. 


























(a) 0 < rij < r„ — no (b) n — no < ni < n + m 

Figure 3: Representation of the distances Vij, Vio in the vulnerability radius 


As described in Section |4.3[ the intensity pcsma of the 
point process representing the transmitters (the MMP) is 
different from the underlying PPP intensity p. Fig. [^depicts 
the value of pcsma in function of p (it corresponds to equation 
( | 11 [ |), we observe that the value of pcsma converges toward 
a maximum value: lim pcsma = ,. 3 „ . It means that, 

p —>-oo 

when the intensity of nodes increases, the intensity of inter- 
ferers converges. In the remainder of the paper, we call this 
situation saturated condition. The convergence implies that 
even in the 2D case where all the nodes are projected on the 
Euclidean plane, the intensity of interferers will not grow as 
much as in the PPP case. 
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Figure 5: Intensity of the MMP in function of the intensity 
of the underlying PPP 


4.6 Comparison with simulation results 

In this section, we compare the 3D theoretical model with 
the 2D theoretical model and with simulation results. The 
goal, as in Section |3.1[ is to evaluate the adequacy of 2D 
and 3D models to represent the probability of coverage for 
different node parameters and building heights. In this sec¬ 
tion, we experiment with two simulation setups: the first 
with 802.11 node parameters, and the second with node pa¬ 
rameters representing 802.15.4 radios. 


The first simulation setup is as follows: the underlying 
PPP is generated in box of 2000m x 2000m x Z with Z € 
20m, 200m, 2000m. We use a larger box than for the PPP 
case because otherwise we do not have enough nodes in the 
MMP after the thinning of the PPP through the contention 
process. The contention process consists in nodes picking a 
random mark uniformly in [0,1]. Then each node checks if it 
can detect neighbors with a lower mark. A node is detected 
if the signal received is above a threshold (the signal strength 
depends on the transmission power, the distance, and on a 
fading random variable). If no detectable neighbor has a 
lower mark, the node is kept in the MMP. 

As in Section [ 3 . 1 1 the 3D intensity of the process is calcu¬ 
lated based on the 2D intensity: p — hi this case, we take 
A = 7 X 10“® so that we have enough nodes in the box and 
the thinning process is not too long (it is notably very long 
for the parameters considered in Section[^. We will see that 
this value is important only if we are not in saturated con¬ 
dition (see the comment on Fig. [^in the previous section). 
For each realization of the MMP, we select the closest node 
to the center of the box to be the emitter, we than compute 
the SIR at no meters from the emitter and check if it is over 
or below the reception threshold. We repeat this process 
1000 times for every d. The parameters of the node and 
the channel are as follows: Pt = lOOmW, Td = —76dBm, 
P = 10, a = 4, p — 1, ed = 10“® and £„ = 10“^. 

Fig. [^depicts the simulation results with Z = 2000m and 
2D (solid curve) and 3D (dashed curved) model predictions 
for the probability of coverage. In this case, the 3D model is 
not in saturated condition. We observe that the 2D and 3D 
models yield very close values. Nevertheless, the simulation 
results closely follows the 3D model predictions. 

Figs. and |7b| present the same comparison as Fig. 
but for more reasonable simulated Z: 200m and 20m. In 
these cases, the 3D model is in saturated condition and we 
observe as in Fig. that the theoretical predictions of the 
2D and 3D models differ. Interestingly, for both Z — 200m 
and Z = 20m the simulation results are closer to the 2D the¬ 
oretical predictions. We note that for Z = 20m the match 
with the 2D model is better than for Z = 200m notably 
for the tail of the curve. These results can be explained 
by the protection radius granted to the emitter thanks to 
the contention process. Indeed, with the considered param- 












(a) Z = 200m (b) Z — 20m 

Figure 7: Comparison of theoretical and simulation results for the MMP model 




Figure 6: Comparison of theoretical and simulation results 
for the MMP model for Z = 2000m 


eters, the average detection range is 251m. Interferers are 
thus very rarely situated above or below the receiver in this 
context (Z = 200m and Z = 20m) and so the vertical dis¬ 
tances do not have a significant impact on the probability 
of coverage. We believe this account for the better match of 
the 2D model. 

The second simulation setup is similar to the first except 
for Td — —GOdBm and Pt = ImW (realistic for a typical 
WSN low power radio [^) and the simulation box size is re¬ 
duced to 200m x 200m x Z with .Z G 20m, 50m, 75m, 100m, 200m. 

Fig-i shows a comparison of the simulation results with 
2D and 3D models for the second simulation setup. We are, 
again, in saturated conditions. In this case, we observe that 
the simulation results are distributed between the 2D and 
3D theoretical models. When Z increases the 3D model is a 
better match for the simulations results, while the 2D is bet¬ 
ter for smaller Z. In this case, the average detection range 
is 32m. From the results presented in Fig. we conclude 
that the 3D model is a good representation of CSMA access 
in dense urban environment if the building height is at least 
two times greater than the average detection radius. 


Figure 8: simu 


In this section, we highlighted the fact that in the CSMA 
case, the model has to be chosen carefully. Indeed, the 3D 
model seems to be relevant only if the average detection 
radius is smaller than the building height which is not the 
case for 802.11. Nevertheless, this model seems satisfying for 
low power low range communications for which the detection 
range is much smaller. 

5. CONCLUSION AND FUTURE WORKS 

In this paper, we investigate the probability of coverage in 
the 3D case for two models: PPP and MMP. We show that 
abusively considering 2D networks when they are 3D can 
lead to either overestimating or underestimating the proba¬ 
bility of coverage depending on the model parameters. We 
also notice interesting differences in the integral forms of the 
expressions when going from 2D to 3D which notably affect 
their tractability. 

By comparing the models with simulations, we observe 
that the 3D model is relevant in the PPP case even if the 
height of the simulation box is small compared to its length 
and width. In the MMP case for the modeling of CSMA, 
it depends on the detection radius. If the detection radius 






























is larger than the building height, the 2D model is closer to 
the simulations, otherwise, the 3D model is more relevant. 

Beside the results on stochastic geometry, this work can 
be seen as an incentive to use 3D model and 3D simulation 
layouts for the design and evaluation of protocol which is 
seldom done in reality. Section [4.6| suggests that it is espe¬ 
cially important for the design of low range loT protocols in 
dense urban areas. 

In the future, we plan to compare the model with traces 
collected in dense urban areas in order to validate our the¬ 
oretical predictions and simulations. Moreover, it would be 
interesting to introduce a fourth dimension in the process in 
order to account for the packet arrivals which are spread in 
time. 
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